Finite-temperature Mott transitions in multi-orbital Hubbard model 



in 
o 
o 

(N 

c 



■4— » 

o3 



C 

o 
o 



> 
o 
in 

O 

in 
o 
-i— > 



T3 

C 

O 

o 



X 



Kensuke Inaba, Akihisa Koga, Sei-ichiro Suga, and Norio Kawakami 
Department of Applied Physics, Osaka University, Suita, Osaka 565-0871, Japan 

(Dated: February 2, 2008) 

We investigate the Mott transitions in the multi-orbital Hubbard model at half-filling by means 
of the self-energy functional approach. The phase diagrams are obtained at finite temperatures for 
the Hubbard model with up to four-fold degenerate bands. We discuss how the first-order Mott 
transition points U c i and U C 2 as well as the critical temperature T c depend on the orbital degeneracy. 
It is elucidated that enhanced orbital fluctuations play a key role to control the Mott transitions in 
the multi-orbital Hubbard model. 

PACS numbers: 71.30.+h, 71.10.Fd 



I. INTRODUCTION 

The Mott transition has been one of the most at- 
tractive topics in strongly correlated electron systems i£ 
There are a number of prototype materials in transi- 
tion metal oxides. A well-known example is V203^£ 
which exhibits the first-order metal-insulator transition 
at a certain critical temperature. The phase diagram 
obtained by systematic experiments is consistent with 
theoretical studies of the single orbital Hubbard model 
by means of dynamical mean field theory (DMFT). 5,6 
This implies that the DMFT treatment of the Hub- 
bard model, which properly takes into account local 
electron correlations, captures the essence of the Mott 
transition. However, in order to give more quantita- 
tive discussions about the Mott transitions, it is indis- 
pensable to incorporate the effect of the orbital degen- 
eracy, which often gives rise to rich phase diagrams, 
as observed for Lai_ x Sr x Mn0 3f £ RTiO^&, etc. More 
recently, the orbital-selective Mott transition has at- 
tracted considerable attention?i2*iiii2ii2ii2ii£ii& in con- 
nection with the materials such as Ca2- x Sr x Ru04i£ and 
La n+1 Ni n 3n +i^±2. 

These interesting experimental findings have stimu- 
lated a number of theoretical works on the Mott tran- 
sitions in the two-orbital Hubbard model by means 
of DMFTiiSiiiiiiiSiiSiSSiSiiSSiSiMiSSiSSiSLSSiS^i^SiSiiS Al- 
though the two-orbital model has been studied in detail, 
a systematic study on the systems having more orbitals 
is still lacking at finite temperatures. For example, the 
finite-temperature Mott transitions for the multi-orbital 
Hubbard models has not been investigated quantitatively 
by DMFT. One of the difficulties lies in the practical cal- 
culation within the DMFT framework. Quantum Monte 
Carlo simulations, which can be a powerful numerical 
method to treat local correlations in DMFT, encounter 
sign problems. Also, Wilson's renormalization group 
method gets more difficult to apply as the number of 
orbitals increases. A much more simplified approach, the 
two-site DMFTjS is not efficient enough to treat finite- 
temperature properties. It is thus desirable to systemati- 
cally investigate electron correlations in the multi-orbital 
Hubbard model at finite temperatures. 

Motivated by the above hot topics, we consider the 



Mott transitions in the multi-orbital Hubbard model at 
zero as well as finite temperatures. For this purpose, 
we make use of a self-energy functional approach (SFA) 
proposed by Potthoff recently^^ This method, which 
is based on the Luttinger-Ward functional approach, 36 
gives a powerful tool to discuss electron correlations. A 
remarkable point is that this approach provides an effi- 
cient way to deal with finite-temperature properties of 
the multi-orbital Hubbard model. The main purpose of 
the present paper is to determine the finite-temperature 
phase diagram of the Hubbard model with up to four- 
fold bands, and quantitatively discuss how the orbital 
degeneracy affects the Mott transitions. 

The paper is organized as follows. In the next section, 
we introduce the model Hamiltonian and briefly summa- 
rize the formulation of SFA. In Sec. IIIII we present the 
detailed analysis based on SFA by exploiting the two- 
orbital model as a prototype system. Then in Sec. IIVI 
we determine the phase diagram, and discuss the nature 
of the Mott transitions in multi-orbital systems in detail. 
A brief summary is given in Sec. ]V\ 



II. MODEL AND METHOD 

We consider a correlated electron system having M 
degenerate orbitals, which is described by the following 
multi-orbital Hubbard Hamiltonian, 



TL = Tin + TL' , 
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where c\ aa (a aa ) creates (annihilates) an electron with 
spin cr(=t, i) and orbital a(= 1, 2, 3 • • • , M) at the i th 
lattice site, t a denotes the hopping integral for orbital a, 
U(U') is the intra-orbital (inter-orbital) Coulomb inter- 
action. In the following, we mainly consider the case of 
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U = U', and give brief discussions for more general cases 
including the Hund coupling in the end of the paper. 

In order to address the Mott transitions in the multi- 
orbital Hubbard model, we use SFA^S We briefly sum- 
marize the essence of SFA. For a correlated electron sys- 
tem, the grand potential is generally expressed as 



fi[E] - F[H] + Trln[-(Go 



(4) 



where F[E] is the Legendre transformation of the 
Luttinger-Ward potential $[G], and G (Go) and E are 
the full (bare) Green's function and the self-energy, re- 
spectively. The condition imposed on the functional 
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E. An impor- 



gives the Dyson equation G~ 
tant point is that the functional form of F[S\ does not 
depend on the detail of the Hamiltonian TLq as far as 
the interaction term TL'(JJ) keeps its shape unchanged. 
This fact allows us to introduce a proper reference sys- 
tem having the same interaction term, which we denote 
as Tircf = Wo(t') + W(U). Then the grand potential is 
written as, 



fi[E(t')] = n(t') 

+ TUnt-CGo^^-SCt'))- 1 ] 
- Trlnr-CGoCtO^-ECt'))- 1 ] 



(6) 



where O(t') and E(t') are the grand potential and the 
self-energy for the reference system, whose bare Green's 
function is denoted as Go(t') -1 = oj + fi — t' (/x: chemical 
potential). The variational condition JSJ is rewritten as 



dt' 



= o. 



(7) 



By choosing the parameters t' for the reference system 
so as to satisfy the condition Q, we can find the best 
system within the reference Hamiltonian, which can ap- 
proximately describe the original correlated system. 

It should be noticed here that SFA provides us 
with an efficient and tractable way to deal with finite- 
temperature properties of the multi-orbital system, 
where standard DMFT combined with numerical tech- 
niques faces difficulties in a practical computation when 
the number of orbitals increases. Another notable point 
in this approach is that the critical behavior can be dis- 
cussed more precisely than the DMFT analysis^ when 
one chooses the same type of the effective impurity 
model. In fact, by comparing the results obtained by 
the two-site DMFTf2&21 the critical point for a single 
orbital Hubbard model obtained by SFA with the two- 
site model 34 is in good agreement with that obtained by 
DMFT with the aid of numerical techniques. 38 ' 39 ' 40 

To discuss the multi-orbital Hubbard model 10 , we ex- 
ploit the Anderson impurity model as a reference system 



in SFAi^ 4 * 3 ^ The Hamiltonian for the reference system is 
explicitly given as 



(8) 
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where ife21(ife2 CT ) creates (annihilates) an electron with a 
spin and a orbital at the k(= 1, 2, • • • Nb) th site, which is 
connected to the i th site in the original lattice. Since we 
consider the Mott transitions without symmetry break- 
ing, we set Eka and Vk a site-independent. Note that in 
the limit of Nb — » oo, Eq. Q is equivalent to the self- 
consistent equation in DMFT. 34 Since the Green's func- 
tion and the self-energy are diagonal with respect to the 
site, spin, and orbital indices, the grand potential per site 
reads 
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+ 2 EE ckjf(u>)6[G' ka (u)}, (10) 

a k=0 

G' 0a (Lu) = [uj + ii-eoa- AaH-E^w)]- 1 , (11) 
G' ka (Lu) = (w + M-ejta) -1 , (k = l,2,.-.,JV 6 ),(12) 

N b 

A a (w) = ^^GLM, (13) 
fc=i 

where R a (cu) = J^p^dz, f(tu) = 1/(1 + e"^) and 
9(u>) is a step function. The grand potential and the self- 
energy for an impurity system are denoted as ili mp 
and S Q (w). 

In the following, we focus on the half- filled Hubbard 
model (/j, = (M — 1/2)U) to discuss how the orbital de- 
generacy (M — 1, 2, 3, 4) affects the Mott transitions. We 
take the simplest reference system (Nb = 1), where we fix 
the parameters (e 0a , E\ a , V\ a ) = (0, fi, V) for each orbital 
a. This simplified treatment is somewhat analogous to 
the two-site DMFT at zero temperature. However, the 
present approach enables us to treat finite-temperature 
quantities systematically, which is a notable advantage 
beyond the two-site DMFT. The condition J7J is now re- 
duced to dQ/dV = 0. We note that the hopping integral 
V between a given site on the original lattice and a ficti- 
tious site corresponds, roughly speaking, to the renor- 
malized bandwidth of quasi-particles. In fact, Fermi- 
liquid properties at zero temperature are determined by 
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the value of V. For instance, small V implies the for- 
mation of heavy quasi-particles, and the system enters 
the insulating phase at V = 0. By calculating the effec- 
tive hybridization V, we thus discuss the stability of the 
metallic state in the multi-orbital systems. 



III. DETAIL OF CALCULATIONS: 
TWO-ORBITAL SYSTEM 

Let us start with the two-orbital Hubbard model, 
by which we illustrate the basic procedure of SFA to 
determine the phase diagram for more general multi- 
orbital cases. We adopt here a semielliptic density of 
states p(e) = (ir/W)y/i - (2e/W) 2 with the bandwidth 
W = 4. 

We first look at the ground-state properties by exam- 
ining the stationary points in the grand potential. In 
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FIG. 1: The grand potential fi(V) for the two-orbital model 
(M — 2) as a function of V for different U at T — 0. Crosses 
correspond to stationary points. 

Fig. ^ we show the grand potential at zero temperature 
for several values of the Coulomb interaction U. When 
U is small, the grand potential has a minimum at finite 
V. Thus the effective bandwidth is finite, stabilizing the 
metallic ground state. It is seen that the stationary point 
moves towards the origin continuously with increasing 
U. For large U, the minimum is located at the origin, 
indicating that the Mott insulating state is stabilized by 
strong electron correlations. In the vicinity of the critical 
point U c , the grand potential can be expanded in powers 
of V, 



Q(V) = ft(0) + AV 2 + 0{V A ). 



(14) 



Therefore, the critical point separating the metallic and 
the insulating phases is characterized by the condition 
A = 0j24 By solving this analytically, we obtain the self- 
consistent equation for the critical point U c as, 

100 f° 

Uc = - 39-y [z-t(z,U e )]p{z)dz 
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FIG. 2: (left panel) The grand potential Q as a function of the 
effective hybridization V for U — 7. Closed circles and trian- 
gles represent the stationary points for the metallic and the 
Mott insulating states. Diamonds denote unstable stationary 
points, (right panel) Stationary values of the grand potential 
as a function of the temperature T for (7 = 7. Symbols are 
the same as in the left panel. 
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■p(z)dz(U) 



with U) = Vc7 2 + z 2 . The detail of the derivation 
is shown in Appendix A. The critical point U c ~ 9.217 
thus obtained is more accurate than the results (U c = 10) 
estimated by the two-site DMFT method4i 

We now consider the competition between the metallic 
and the Mott insulating phases at finite temperatures, 
following the way for the single band case^ The grand 
potential is shown in Fig. [2 at finite temperatures. It is 
seen that two minima appear at low temperatures. One 
of the minima is located at larger V , which corresponds to 
the metallic solution, since it is continuously connected 
to the metallic one at zero temperature. The other is 
adiabatically connected to V = at zero temperature, 
so that this solution characterizes the Mott insulating 
state. Such a double-well structure at low temperatures 
causes the first-order transition accompanied by hystere- 
sis. In fact, as increasing temperature, the stationary 
point for the metallic state disappears around T c i = 0.05, 
where the Mott transition occurs to the insulating phase. 
On the other hand, as decreasing temperature, the Mott 
insulating phase realized at high temperatures is stable 
except for zero temperature, since the corresponding lo- 
cal minimum always exists at finite temperatures. The 
first-order transition temperature T c = 0.045 for (7 = 7 
is determined by the crossing point of the two minima in 
the grand potential, as shown in Fig. [5] This hystere- 
sis induced by the first-order transition is also observed 
when the interaction U is varied. We show the effective 
hybridization V at T = 0.04 in Fig. |3| Starting from 
the metallic state, the increase of the Coulomb interac- 
tion decreases the effective hybridization V, and finally 
triggers the Mott transition to the insulating phase at 
U C 2 = 7.2, where we observe the discontinuity in V . On 
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FIG. 3: The effective hybridization V as a function of U at FIG. 4: The quasi-particle weight Z as a function of the 

T = 0.04. Inset shows the grand potential as a function of Coulomb interaction U for the system with M = 1, 2, 3,4 at 

17 at T = 0.04, where fi M (fir) is the grand potential for the T = 0. The inset is the enlarged figure in the small U region, 
metallic (insulating) state. 



the contrary, the Mott transition occurs at U c \ = 6.7 
when the interaction decreases. In this parameter regime, 
the level crossing in the minima of the grand potential 
appears at U c = 7.1 (inset of Fig. EJ), which defines the 
first-order transition point at T = 0.04. 

We have seen here that the two-orbital Hubbard model 
exhibits qualitatively similar properties to the single- 
band cas o 5 ! 34 as far as the nature of the Mott transi- 
tions is concerned. In the following section, we give more 
quantitative discussions about the Mott transitions by 
using the phase diagrams obtained for the multi-orbital 
Hubbard model with M = 1 ~ 4. 



IV. PHASE DIAGRAMS OF MULTI-ORBITAL 

SYSTEMS 

We now determine the phase diagrams of the multi- 
orbital Hubbard model. We first calculate the renormal- 
ization factor Z = (1 — dT,(aj)/duj)~ 1 , which is inversely 
proportional to the effective mass, to characterize the 
metallic ground state at zero temperature. The results 
are shown in Fig. The introduction of the Coulomb 
interaction decreases the renormalization factor, making 
the mass of quasi-particles heavier. In the small U region, 
there is not a big difference in the behavior of Z among 
different M cases. If we look at Z in more detail (inset of 
Fig. 0J, the mass gets slightly heavier (smaller Z) as M 
increases. This comes from the fact that the electron cor- 
relations are somewhat enhanced for multi-orbital cases 
by the inter-orbital Coulomb interaction. On the other 
hand, further increase of the Coulomb interaction leads 
to quite different behavior: as the number of orbitals 
M increases, the mass is less renormalized, which makes 
the metallic state more stable up to large U 4242iii This 
implies that the inter-orbital interaction enhances orbital 



fluctuations for large U, which in turn stabilize the metal- 
lic state, as pointed out in the two-orbital case4i In fact, 
the critical points U c2 (T = 0) = 9.2173,12.6044 and 
15.9958 for M — 2,3 and 4 (see Appendix), are much 
larger than U c = 5.84 for the single-orbital models 
Therefore, enhanced orbital fluctuations play a key role 
to stabilize the metallic state at zero temperature. 

We now move to the finite-temperature properties. 
The phase diagrams obtained in the way mentioned in 
the previous section are shown in Fig. There are three 
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FIG. 5: Phase diagrams for the degenerate Hubbard model 
with the orbital degeneracy M = 1,2, 3, 4. 

phase boundaries for each system, U&, U c i and U c . The 
double-well structure in the grand potential gives rise 
to two kinds of transitions with the discontinuity in the 
physical quantities. On the phase boundary U c i(U C 2), 
the Mott transition occurs from the insulating (metallic) 
state to the metallic (insulating) state as U decreases 
(increases). The area surrounded by U c i and U C 2 is re- 
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ferred to as the coexistence phaseA24 The critical point 
U c is determined so that the two minima of the grand 
potential take the same value. At zero temperature, the 
critical points U c and £/ c2 merge to give the continuous 
Mott transition, because the double-well structure dis- 
appears at T = 0. Therefore, the introduction of the 
temperature drastically changes the phase boundary as 
discussed in the single-orbital caseA^ It is found that 
as temperature increases, the phase boundaries U c %, U C 2, 
and U c merge at the critical temperature T c , where the 
second-order phase transition occurs. 

The above characteristic properties in the Mott transi- 
tions are qualitatively the same as the single-orbital case. 
We now discuss the effects due to orbital fluctuations 
quantitatively. First note that the coexistence region 
bounded by the phase boundaries U& and U C 2 is enlarged 
when the number of orbital M increases. Furthermore, 
we see that the M-dependence of the critical points U& , 
U C (T C ) and the critical temperature T c is different from 
that of U C 2- To clarify this point, the critical values are 
plotted as a function of M in Fig. It is found that 
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FIG. 6: The critical points U C (T = 0),U C (T = T B ), U c i(T -» 
0) and the critical temperature T c as a function of the orbital 
degeneracy M. 

U C 2 is proportional to the orbital degeneracy M, while 
Ud and U C (T C ) the square root of M. These results are 
consistent with the previous results However, 
we wish to note that the hnite-temperature properties 
such as the phase diagrams have not been clarified quan- 
titatively so far. In particular, the results on the in- 
dependence of the critical temperature T c obtained here 
are beyond the qualitative discussions by Florens et al.^* 
who obtained only the upper bound for T c cx M . 

To conclude this section, we would like to briefly dis- 
cuss the effects of the Hund exchange coupling. For this 
purpose, we add the following term to the Hamiltonian 
©, 

i a.<a' tr 

- J Y1 ( c Lt c «*i c L'i c «*'t + H - c -) 

i a<a' 

~ J Y1Y1 ( C !aT C Ial C »'T c M 'i + H.C.) (16) 
i a<a' 




FIG. 7: Phase diagrams for the two-orbital Hubbard model 
with the Hund coupling J = 0.1(7 under the condition U = 
U' + 2 J. For clarity, we also show the results at J = 0. 

with J > 0, where we impose the condition U = U' + 2J 
due to the symmetry requirement. Since the calculation 
gets somewhat difficult in the presence of J, we show the 
results obtained for the representative case of M = 2. 
The obtained phase diagram is shown in Fig. and 
compared with the J =0 case. First, we notice that 
upon introducing J, the phase boundaries are immedi- 
ately shifted to the weak-interaction regime, and there- 
fore the metallic state gets unstable for large U. Corre- 
spondingly, the coexistence region surrounded by the first 
order transitions shrinks as J increases. This remarkable 
tendency indeed reflects the fact that the metallic state is 
stabilized by enhanced orbital fluctuations in the case of 
U = U': the Hund coupling suppresses such orbital fluc- 
tuations, and stabilizes the Mott-insulating phase. An- 
other point to be mentioned is that the Mott transition 
becomes of first-order even at zero temperature in the 
presence of J^^l^ since the subtle balance realized at 
T = in the case of U = U' is not kept anymore for 
finite J. Since there is another claim that the second 
order transition is possible for certain choices of param- 
eters at T = 0j3a which may depend on the strength of 
Hund coupling, we need further detailed discussions to 
clarify this problem. Although we have presented only 
the results for the M = 2 case, we expect that the effects 
of J on the phase diagram should be essentially the same 
as shown here: the introduction of J suppresses orbital 
fluctuations, thus favors the Mott insulating phase even 
in the small U regime. It also shrinks the coexistence 
phase dramatically. 



V. SUMMARY 

We have investigated the Mott transitions in the multi- 
orbital Hubbard model at zero and finite temperatures by 
means of the self-energy functional approach. By choos- 
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ing the Anderson impurity model as a reference system, 
we have discussed how the orbital degeneracy affects the 
nature of the Mott transitions. We have obtained the 
finite-temperature phase diagram for the system having 
M(< 4) degenerate bands. Although the phase diagrams 
show qualitatively similar features irrespective of the or- 
bital degeneracy, there are some remarkable effects due to 
degenerate orbitals. In particular, it has been shown that 
enhanced orbital fluctuations make the metallic phase 
more stable even at finite temperatures. Therefore, if 
such fluctuations are suppressed, the metallic state is ex- 
pected to be unstable, which has been shown to be the 
case by considering the system with the Hund coupling. 
Also, it has quantitatively been clarified how distinctly 
the M-dependence appears in the critical points U c i and 
U C 2'. the critical point U c \ depends on the square root 
of the orbital degeneracy M, while the critical point U C 2 
is proportional to M. This analysis concludes that the 
critical temperature T c is proportional to the square root 
of M, which may be important to understand the Mott 
transitions in the real materials with orbital degeneracy. 

The self-energy functional approach used in the paper 
allows quantitatively reliable discussions for multi-orbital 
systems at finite temperatures. Since this formalism pro- 
vides a tractable way to incorporate spin and charge or- 
dered states, which have been neglected in this paper, 
it may be used for more detailed study of the finite- 
temperature properties in multi-orbital Mott systems. 
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APPENDIX A 

We analytically determine the critical value of U for 
the Mott transition in the multi-orbital Hubbard model 
at zero temperature. The grand potential Eq. (|10fl is 
rewritten as 

n/L = n tap + 2MJ2 J dzp(z)wi(z)0[-Wi(z)] 

i 

- 2M^> i 0(-u;0, (Al) 

i 

where f2i mp is the grand potential of the reference sys- 
tem and 0(x) is the step function. uj[,uji{z) are the poles 
of the impurity Green's function of the reference system, 
G'(v), and the approximated Green function of the orig- 
inal system, G(w; z) = 1 / (u> + p — z — H(ui)) , respectively. 



First we consider the two-orbital system. In the atomic 
limit V = 0, the Green's function and the self-energy at 
the impurity site for the reference system are given by 



Gatom( w ) 
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lu - U/2 
IP 
' 4w' 



U/2 



(A2) 
(A3) 



where the poles of G' &tom are ±U/2. Then the Green's 
function of the original lattice model is given by 



U 2 

. M , m (u);z) = UJ-Z+-, 



G- 1 



which has the poles at 

u±(z) = - (z ± V~z 



U 2 



(A4) 



(A5) 



As discussed in the text, we can expand the quantities 
in V around the atomic limit. The grand potential is 
il' = -2U - 2W 2 /U + 0{V 4 ). We also analyze the 
Green's functions G'(u>) and G{u)\z) around the atomic 
limit and obtain their poles up to the second order in 
V. Here we need only the poles in the negative energy 
region, because the poles in the positive energy region do 
not affect fi, see Eq. (|A1|) . The Green's functions G\ui) 
and G(oj; z) have eight poles in the M = 2 case, and thus 
four poles in the negative region. The negative poles of 
G'(lj) are 
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and the negative poles of G{oj\ z) are 

wi(z) = W0zV 2 /U 2 + (D(V 4 ), 

lu 2 (z) = lu-{z) + B(z/U)V 2 /U + 0(V 4 ), 

lu 3 (z) = -U/2-30V 2 /U + 0(V 4 ), 



cu 4 (z) = —217 — 



168z + 390C7 
15U + 8z 



V 2 /U + 0(V 4 ), 



(A6) 



where 
B(x) = 



6525 + 8714x 2 - 320Cte 4 50x(237 - 64x 2 ) 



(225 - 64x 2 )V 1 + x 2 



225 - 64x 2 



(A7) 



The grand potential is written down as 

/•CO 

fl/L = -4 / dzp{z)S,(z 7 U) 
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dzp(z) (z-Z(z,U)) 



dzp(z) 



Zl7U + 83£(z,U) 
[17U + 8£(z,U)]S(z,U) 



The 



-o(v 4 ), 

(A8) 
condition 



where U) 
d 2 fl/dV 2 \v=o — is satisfied at the Mott transi- 
tion point U c m^ Given that the free density of state is 
symmetric p(z) = p(—z), we derive the self-consistent 
equation for U c , Eq. (|15|l . 

Similarly, we can analyze the cases of M — 3 and M = 
4. The self-consistent equation for U c at M — 3 reads 



49 



U c = / [z-t(z,U D ))]p(z)dz 



_Uc 

80 

and at M = 4 is 

L/ c = 

135 
135 



545C/ C + 80£(z, C/ c ) 

[nu c + 8£(z, u c )]£(z, u c 



-p(z)dz, 



(A9) 



[z - f («, C/ c )] p(«)dz 

837t/ e + 63g(z, E/ c ) 
^ [17E/ C + 8£(z, U c )]£(z, U c ) 



p(z)dz. 
(A10) 



The solution of these equations gives the accurate values 
of the critical point: U c = 12.6044 and 15.9958 for M = 3 
and 4, respectively. 
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